Gravitational signature of Schwarzschild black holes 
in dynamical Chern-Simons gravity 



C. Molina,i'B Paolo Pani,2^[t| Vitor Cardoso, 3- and Leonardo Gualtieri^'d 

' Escola de Artes, Ciencias e Humanidades, Universidade de Sao Paulo, 
Avenida Arlmdo Bettio 1000, 03828-000 Sdo Paulo, SP, Brazil 
Dipartimento di Fisica, Universita di Cagliari, and INFN sezione di Cagliari, Cittadella Universitaria 09042 Monserrato, Italy 

Centro Multidisciplinar de Astrofisica - CENTRA, Dept. de Fisica, 
Instituto Superior Tecnico, Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal 
"^Department of Physics and Astronomy, The University of Mississippi, University, MS 38677-1848, USA 
^Dipartimento di Fisica, Universita di Roma "Sapienza" and Sezione INFN Romal, P. A. Moro 5, 00185, Roma, Italy 

Dynamical Chern-Simons gravity is an extension of General Relativity in which the gravitational 
field is coupled to a scalar field through a parity- violating Chern-Simons term. In this framework, we 
study perturbations of spherically symmetric black hole spacetimes, assuming that the background 
scalar field vanishes. Our results suggest that these spacetimes are stable, and small perturbations 
die away as a ringdown. However, in contrast to standard General Relativity, the gravitational 
waveforms are also driven by the scalar field. Thus, the gravitational oscillation modes of black 
holes carry imprints of the coupling to the scalar field. This is a smoking gun for Chern-Simons 
theory and could be tested with gravitational-wave detectors, such as LIGO or LISA. For negative 
values of the coupling constant, ghosts are known to arise, and we explicitly verify their appearance 
numerically. Our results are validated using both time evolution and frequency domain methods. 
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I. INTRODUCTION 

In Chern-Simons gravity [IHll the Einstcin-Hilbert ac- 
tion is modified by adding a parity-violating Chern- 
Simons term, which couples to gravity via a scalar field. 
This correction could explain several problems of cosmol- 
ogy [3i8|. Furthermore, a Chern-Simons term arises in 
many versions of string theory and of loop quantum 
gravity [Tol - [l^ . and Chern-Simons gravity can be recov- 
ered by truncation of low energy effective string models 

EH. 

When Chern-Simons gravity was first formulated, the 
scalar field was considered as a prescribed function. Later 
on, it was understood that this a priori choice is not re- 
ally motivated (see the discussion in Ref. [l^)- Then, 
dynamical Chern-Simons (DCS) gravity has been formu- 
lated , in which the scalar field is treated as a dynam- 
ical field. 

Since DCS gravity has a characteristic signature (the 
Chern-Simons term violates parity), there is the exciting 
prospect of testing its predictions against astrophysical 
observations. This has motivated a large body of work 
on the subject (for a review on DCS gravity and its astro- 
physical consequences see Ref. [iBl)- In this context, the 
study of black hole (BH) perturbations is very promis- 
ing, since astrophysical black holes are probably the most 
appropriate objects to probe the strong field regime of 
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General Relativity ■ 

The first study of BH perturbations in the context of 
DCS gravity has been carried out in Ref. [l^, where it 
was found that, if the background solution contains a 
(spherically symmetric) scalar field, polar and axial per- 
turbations of DCS BHs are coupled, and the equations 
describing them are extremely involved. Recently, in 
Ref. [l^ (hereafter. Paper I), some of us found that, when 
the background scalar field vanishes, polar and axial 
gravitational perturbations of a Schwarzschild BH decou- 
ple, and only axial parity perturbations are affected by 
the Chern-Simons scalar field. We also found that under 
this assumption the gravitational and scalar perturba- 
tions are described by a coupled system of two second or- 
der ordinary differential equations (ODEs). The numer- 
ical integration of this system to find the quasi-normal 
modes (QNMs) of Schwarzschild DCS BHs is challenging, 
due to the same asymptotic divergence which prevented 
for many years the numerical computation of QNMs for 
Schwarzschild BHs [lO-ill. Therefore, in Paper I the 
QNMs of Schwarzschild DCS BHs were not investigated 
thoroughly. It is remarkable that there are very few 
studies of this kind of system, i.e., QNMs described by 
coupled ODEs (one interesting work is presented in Ref. 
[2^). In Paper I we also claimed that Schwarzschild DCS 
BHs are unstable for a specific range of the parameters 
of the theory. This result was the consequence of a sign 
error in the derivation of the perturbation equations; on 
the contrary, as we discuss in this paper, there is strong 
evidence that these spacetimes are stable. 

In this paper we complete the study of Schwarzschild 
DCS perturbations, performing a thorough numerical 
analysis of the perturbation equations. We employ two 
different - and completely independent - numerical ap- 
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preaches: time evolution and a formulation of the fre- 
quency domain approach (25| which has never been ap- 
plied before to the study of instability in black hole space- 
times. The results of the two independent methods agree 
very well, typically within an accuracy of 0.1%, validat- 
ing each other. 

The main result we find is that any perturbation de- 
cays at late-time as a damped sinusoid. This is known as 
the ringdown phase, where the black hole radiates all ex- 
cess hairs in its lowest QNMs [22I, [2^ . What is new here, 
and with important implications for tests of DCS grav- 
ity, is that the gravitational sector has two distinct sets of 
QNMs. For large values of the constant (3 (associated to 
the dynamical coupling of the scalar field) , these two sets 
coincide with the usual gravitational QNMs and scalar 
field QNMs of General Relativity. This result enables 
simple, yet fundamental tests on DCS gravity. By mea- 
suring (or not) these two different modes, one could effec- 
tively constrain DCS gravity through gravitational- wave 
observations. For instance, detection of ringdown modes 
with a signal-to-noise ratio > 6 (feasible with both the 
Earth-based LIGO and the space-based detector LISA), 
could allow one to test DCS gravity if the mass of the 
BH is known, for instance through observations of the 
inspiral phase of black hole binaries. For signal-to-noise 
ratios > 150 one could be able to discriminate between 
DCS gravity and standard General Relativity without 
any further knowledge of the BH parameters. 

A summary of our results 

For the reader wishing to skip the technical details of 
the rest of the paper, the following is a brief summary of 
our results. 

(i) Two complementary numerical methods were de- 
veloped and employed. They are completely inde- 
pendent and their concordance is very good. 

(ii) For small values of the coupling constant {M*f3 < 
0.5), the perturbative dynamics is characterized by 
a stable exponentially decaying phase. The inter- 
mediate late time evolution is dominated by 



*(i, ''fixed) = 




with Re[wno] = and Im[cJno] < (with our sign 
conventions, a QNM is stable if Im[a;] < 0). Our 
results for the non-oscillatory frequency values arc 
compatible with the expression: 

..„ = -0.04024,MW.«f(l + H:5p5-if?)^ 

(1.2) 

(iii) For intermediate values of M"*/?, field evolution is 
dominated by a stable oscillatory phase. We have 



detected two oscillatory modes, named here "grav- 
itational" and "scalar" modes. Although the time 
profiles of the gravitational perturbation ^I* and of 
the scalar field & are usually different, they consist 
on different superpositions of the same modes. 

(iv) In the /3 00 limit, these "gravitational" and 
"scalar" branches coincide with actual gravita- 
tional and scalar modes of Schwarzschild BHs in 
General Relativity. In this regime, we report that 
for ^ = 2, we find that the gravitational perturba- 
tion oscillates with a combination of the two modes 

Ma;grav(1') = 0.3736-1 0.08899, (1.3) 
MuscW = 0.4837-1 0.09671. (1.4) 

These numbers correspond to the lowest mode of 
pure gravitational and scalar quasi-normal frequen- 
cies in Einstein's theory The scalar field per- 
turbation, instead, oscillates with the mode Wgc 
only. This behavior can be easily understood by 
looking at the form of the equations in this limit. 

(v) At late times, the field decays with a power-law 
tail, as t~(^^+'^'. The tails do not depend on 
/3 or M. Note that the same behavior charac- 
terizes Schwarzschild BHs [26], implying that a 
gravitational- wave observation of the tail would not 
be able to discriminate DCS gravity from General 
Relativity. 

(vi) An extensive investigation of BH oscillations, per- 
formed using two different numerical approaches, 
only yields stable modes, either oscillating or non- 
oscillating. This gives strong indications that 
Schwarzschild BHs in DCS modified gravity are sta- 
ble against axial and polar perturbations. 

(vii) We also discuss how the inclusion of a non- 
vanishing scalar potential in the Lagrangian affects 
the QNM spectrum. We focus on potentials of the 
form 

V{d) = m^d^ + 0{d^) (1.5) 

and find that in the /? — > 00 limit this inclusion 
only affects the scalar branch of QNMs, while the 
gravitational branch is unaltered. When < 
100, also the gravitational sector is affected by the 
scalar potential. 

The paper is organized as follows. In Section [H] we 
briefly review the derivation of the perturbation equa- 
tions in DCS gravity. In Section iHll we describe the time 
domain and frequency domain numerical approaches that 
we have employed to solve the perturbation equations. 
In Section [IV] we present our results in the time and fre- 
quency domains. In Section |Vl a possible observational 
signature of DCS gravity is discussed. Implications and 
final remarks are presented in Section [VII 

In Appendix |^ we discuss ghost-like instabilities aris- 
ing when the wrong sign of the kinetic term in the action 
is chosen, i.e. when /3 < in Eg. l2.1l below. 
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II. PERTURBATION EQUATIONS AND 
DYNAMICAL STABILITY 



The action of DCS gravity is [1 
S ~ K I 



gR+j I d^x^/^^ *RR 



J [g'^^'Va^Vbd + + 5„,at (2.1) 



where -d is the scalar field and 



t TDTD D ^baef T>cd 

KK — -Kabcd^ n , 



(2.2) 



Wc use geometrical units c = G = 1 so that k = I/IGtt. 
Furthermore, we neglect T^(i?) (this assumption will be 
relaxed in Section |IVE|) . and consider the vacuum solu- 
tions (^mat = 0). The equations of motion are 



Rab 



-167raCa6 + 87r/3z9,ai?,& 



ad = *RR 

4/3 



where 



(2.3) 
(2.4) 

(2.5) 



In a spherically symmetric background, *RR = = C"'' 
and Eqs. (|2.3p . (|2.4p reduce to usual Einstein gravity 
minimally coupled to a scalar field 



No-hair theorems (27[ state that the Schwarzschild solu- 
tion, with vanishing scalar field, is the only static spheri- 
cally symmetric solution of the equations above. We then 
consider perturbations of a Schwarzschild BH with a van- 
ishing background scalar field. We expand the gravita- 
tional perturbations in tensor spherical harmonics, build- 
ing the Zerilli and Regge- Wheeler functions. The scalar 
field is expanded in scalar spherical harmonics as 



Eq. (|2.3p implies (see Paper I) that polar parity gravi- 
tational perturbations (described by the Zerilli function) 
are not affected by the Chern-Simons scalar, and then 
the corresponding QNMs are the well-known modes of 
Schwarzschild BHs. Axial parity gravitational perturba- 
tions = iQ^"'/i^ (where Q^™ is the Regge- Wheeler 
function, defined as in Paper I) are instead coupled with 
the scalar field. From here onwards, we will drop the ^™ 
superscripts. 

Eqs. (j2.3p . (|2.4p reduce to the following set of coupled 
ordinary differential equations for the perturbations Q(r) 
and 'l'(r), in terms of which one can completely charac- 
terize the axial parity metric perturbations and the scalar 
(2.6) field respectively: 



^ + 

drl 



f 



mTTAPa^ 



6M 
2M 



^1 



e = 



967rM/ 



(f + 2)!6A/a 
■' {1-2)1 r^p 



(2.8) 
(2.9) 



with /(r) = 1 - 2M/r and = r + 2A/ln(r/2A/ - 1). 
Note that third time-derivatives (i.e. terms proportional 
to w'^) do not arise in the perturbation equations above 
(although they are generally expected from Eqs. (|2.3p - 
p.4p ) because of the vanishing of the background Ricci 
tensor in Eq. <\2.h\ . Therefore the Schwarzschild back- 
ground does not suffer from problems related to ill- 
posedness of the theory, the so-called Ostrogradski in- 
stability (see Refs. [l6l . [28l]). We also remark that the 
instability found in Paper I for /3M^ < 27r was an arti- 
fact of a wrong sign in the definition of *RR, that has 
yield a change in the sign of (3 in the perturbation equa- 
tions. This is equivalent to consider the equations of the 
DCS theory with /3 < 0, which is indeed expected to be 
unstable, as discussed in Appendix Rl 



Re-scaling and the General Relativity limit 

Under the replacement /3 — >■ a^/? and Q — >■ B/a, one 
can set a = 1 in the perturbation equations (|2.8|) and 
p.9|) . which we will hereafter assume. Indeed, as dis- 
cussed in [l6j . the parameters of the theory are redun- 
dant, and it is always possible to fix one of them. 

We remark, however, that there is a subtle formal dif- 
ference between the theory with a, /3 and the theory with 
a — \. Indeed, the General Relativity limit of the former 
is obtained by taking /3 — > oo and a — > 0; the General 
Relativity limit of the latter is obtained by taking /3 — >■ oo 
and by considering the solutions with = 0. In other 
words, once we fix a = 1, General Relativity is not sim- 
ply a limit of the DCS theory: it is a particular subset of 
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the solution space of the /3 — )• od Umit of the theory. 

III. NUMERICAL APPROACH 
A. Time domain evolution 



The system p.Sp . ()2.9p can be written as 

$ = y $ 



dt^ drl 
where we have defined 



$ = 



V = 



(3.1) 



(3.2) 



V21 V22 

and the elements of the matrix potential V are given by 





= fir) 


V12 


= fir) 


V21 


= fir) 


V22 


- fir) 



£{£ + !) 



6M 



967rM 



6M{£ + 2)1 1 
/3(£-2)! ^ 

?+l) 5767rA/2\ 2 A/ 



(3.3) 
(3.4) 
(3.5) 

(3.6) 



Using the light-cone variables u ~ r^, — t and v ~ + 1 
one can write 

A discretized version of Eq. (|3.7p is 

$(Af) - $(i;) - $(VF) + $(S') 
AuAv 

= ^^ViS)ME) + ^W)] , (3.8) 

where the points N, E, W, S are defined as follows: = 
{u + A,v + A), W = {u + A, v), E = {u,v + A) and 
5 = iu,v). With the expression p.Sp . the region of in- 
terest in the u — v plane is covered, using the value of the 
field at three points in order to calculate the fourth one. 
As the integration proceeds, the values of ^(t, rfixod) are 
extracted [ll,!!^]. 

The initial data consist of the expressions on the sub- 
manifolds {u > 0,v = 0) and (u = 0, i; > 0) for the 
vector 



(f>(M, v) 



*(m,u) 
e{u,v) 



(3.9) 



For most of the numerical evolutions presented here the 
initial data have the form 



$(0,i;) 



-(u-Uc) /2cr 



(3.10) 
(3.11) 



with Vc = 10.0 and a = 1.0. 

From results on BH oscillations in General Relativ- 
ity [2l| we expect that the main characteristics of the 
time-evolution profiles (after a transient initial regime) 
are insensitive to the choice of the initial data, provided 
that they are localized. To check if this actually occurs 
in the present case, and rule out any eventual influence 
of initial data on late time results, we have considered 
different choices for the initial data: 



• Gaussian initial data 

g-(lI-«el)^/2(Ti 



$(0,w 



A 



2 e 



^(l'-«c2) /2(T2 



$(w,0) 



(3.12) 

The initial u-functions are localized, with differ- 
ent peaks for the ^ and O components. Although 
strictly speaking they do not have compact sup- 
ports, they are (numerically) zero far away from 
the peaks. 



• Compact support pulses 



mv) = fiv) 



$(u,0) 



(3.13) 



We have chosen two different functions / 
/i(^')i f^iv)- The first choice corresponds to 



fiiv) = 



{v - V2) iv - Vi) 



{V2 - Vi) 



Vi < V < V2 



(3.14) 



and zero elsewhere. This is a localized and smooth 
pulse with a compact support. Our second choice 
corresponds to 



/2(w) = 1 , Vi <V <V2, 



(3.15) 



and zero elsewhere. It is a localized but not con- 
tinuous pulse with a compact support. 

We have verified that the numerical results (after a tran- 
sient regime) do not depend on the initial data. 



B. Iteration scheme in the frequency domain 

We now present an alternative, and complementary, 
numerical method, which is an application of Newton's 
iteration scheme to the shooting method [2^ . 

Let us define ujo as the trial eigenfrequency of the eigen- 
value problem defined by Eqs. (|2.8p . (|2.9p . The corre- 
sponding solutions and 9o satisfy the following set of 
equations 



^'^{r,) + {uj^ -Vn)'foir.) 
e'^ir,) + iu' -V22)eoir.) 



VuSoir,), (3.16) 
V2iMr*), (3.17) 



5 



and it is hereafter understood that all these quantities are 
evaluated at the trial frequency ojq. In order to compute 
QNMs we require the following boundary conditions 



$o(±oo) 



^'o(±oo) 
eo(±oo) 



A± 
B± 



(3.18) 



When Im[w] < Eq. ([3l8| defines (stable) QNMs, while 
when Im[a;] > we have "bound-state-like" boundary 
conditions, i.e. $o ^ at ±03 (see Paper I) and 

the corresponding modes are unstable. The numerical 
method described in the rest of this section is capable to 
find both stable and unstable modes. 

The idea is to "shoot" from each of the boundaries 
to a matching point where the wave functions and their 
derivatives are required to be continuous. In general, ujq 
is not the true eigenfrcquency, and one of the continuity 
equations for and 89 is not satisfied. Without loss 
of generality, we can choose either Qq or Oq to be the 
function which does not satisfy the continuity condition. 
Moreover we consider the matching point to be at = 0. 
Namely we assume that 

M = [[*o]] = [[eo]] = o, [[e^,]]^o, (3.19) 

where we define [[. . . ]] as the difference between the 
limits of the corresponding quantity as — >■ 0±. We 
checked that our numerical results do not depend on 
the choice of the matching point within a wide range 
around = 0. We perform two integrations: one start- 
ing at +00 (numerically, at = ri^^ ^ M) inward to 
r^, — 0, and the other one starting at — c» (numerically, at 



.(2) 



^ ~M) outward to r^, ~ 0. At both infinities, 



we expand solution in series as follows 



Ah 
Bh 

Aoo 




2M)" 



(3.20) 

at r^, <C —M and at ^ M, respectively. In comput- 
ing QNMs this way it is important to choose appropriate 
values of numerical infinities, because numerical instabil- 
ities may arise by considering too large values for ri^'' and 
ri^'' [iOl ■ In fact at both infinities the general solution will 
be a mixture of exponentially growing and exponentially 
suppressed modes and (in order to compute QNMs) we 
must select pure exponentially growing modes. Problems 
arise when too large values for rl^-* and rP'' are chosen, 
because in that case contributions from unwanted expo- 
nentially suppressed modes can be significant after the 
integration, due to numerical errors. This problem can 
be circumvented by choosing small enough values of nu- 
merical infinities, say jri*'! ^ lOM, and by considering 
large enough order of series expansion N, say TV > 10. 
In this way, though ^ 10 A/ is not very large (typically. 



for the modes we find, |10Ma;| ~ 3 — 5) the series well 
approximates the correct solution. This problem does 
not arise in the computation of unstable modes (see Ap- 
pendix |X| , since in that case we simply impose Dirichlet 
conditions at both infinities. 

In order to obtain solutions satisfying continuity condi- 
tions (j3.19p at the matching point, we compute two linear 
independent solutions and wc construct an appropriate 
linear combination of them, which satisfies the required 
conditions. The first solution is obtained by choosing 
A = 1 and a generic value B = Bq in the series expan- 
sion, whereas the second solution is obtained by choosing 
B = 1 and a generic value oi A = Aq. We shall denote 
the first solution as ^+ and the second one as In 
order to have two linear independent solutions we also 
require ^0^0 1- 

The procedure outlined above is adopted twice: once 
for and i?oo and once for Ah and Bh- Accord- 
ingly, we perform four numerical integrations: two from 
ri^^ and two from ri^^ up to = and we obtain 
{^i^^\r^),ei^^\r,)) and {^^g\r^),ef\r^)) respec- 
tively. Finally we construct a linear combination of solu- 
tions: 



60 



= < 



flight 

QHght 
Qloft 



+ b 



right 



right 



eicft 



> 0, 



and we choose a, 6, c, d in order to satisfy the continu- 
ity conditions, Eqs. (|3.19p . The net result of this pro- 
cedure is a set of solutions {\I'o(7v), Bo(?'*)} which have 
the correct boundary conditions and which are continu- 
ous everywhere with ^'q also continuous everywhere. The 
discontinuity [[Qq]] 7^ is related to the choice of a trial 
eigenvalue wq, which is not the correct eigenfrcquency. 

Let us now denote with uji the correction to the trial 
eigenvalue, i.e. lj = ojq + uji. If cji is a small correction, 
i.e. cji <C wq, then [2^ 



UJi 



Jdr^ 



(3.21) 

where we have defined P(r*) = — cj^ + Vu, Q{ri,) ~ 
-u? + Rir*) = V12, S'(r*) = V21. In our case 
dP/dujo = dQ/dojQ = —2ujo and dR/duJo = dS/dujQ = 0. 
Moreover in the Eq. p.2ip Aq and yUo are the solutions of 
the conjugate equations of Eqs. (|3.16p - (|3.17p 



A;,'(r.) + (c.^ 
^o(r^) + {u? 



Vi2^0 ■ 



(3.22) 
(3.23) 



The correction p.2ip has been computed in Ref. |25| for 
the case of "bound-state like" boundary conditions. In- 
terestingly enough, it is also valid for the more general 
case of boundary conditions defined in Eq. p.lSp . In 
fact it is straightforward to show that contributions to 
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Eq. (|3.21[) arising from boundary conditions p.l8|) cancel 
each others out, if the same boundary conditions are also 
imposed on Aq and /iQ. Therefore Eq. p.2ip can be used 
in an iteration scheme until we reach the required accu- 
racy. We find that convergence usually occurs, within the 
required precision (typically |[[6o]] /6o(+0)| < 10"'^), in 
less then 50 iterations. However, we cannot find the 
entire QNM spectrum using this method. Indeed even 
the single equation version of this method fails to find 
first overtones of Schwarzschild BHs in General Relativ- 
ity . This is the reason why, as discussed in the next 
section, we can find QNMs with this approach only for 
M^f3 > 0.5. For smaller values of /3 the iteration scheme 
ceases to converge. 




2 3 4 5 6 7 8 9 10 11 12 

; 



IV. NUMERICAL RESULTS 

In this Section we present the results of our numeri- 
cal integrations, performed using both the time domain 
approach and the iteration scheme approach in the fre- 
quency domain. The results for time domain evolu- 
tions refer to Gaussian initial data, with a Gaussian 
wave-packet characterized by Vc = 10.0 and cr = 1.0 in 
Eq. dSH]); the field is extracted at = 50.0A/. 

A. Small M*/3 limit 

For small values of M^f3 (< 0.5), the perturbative dy- 
namics is characterized by a stable exponential mode 
phase. The intermediate late time evolution is dominated 

by 



*(i,?'fixcd) = e"^ 



(4.1) 



with Re[a;„o] = and Im[ti;,io] < 0. 

After an extensive numerical exploration performed us- 
ing the time domain approach, the non-oscillatory fre- 
quencies Wno obtained are consistent with the expression 



Mojno = -0.04024(M''/3)°-44^ 
which is illustrated in Fig. [T] 



1-f- 



2.0953 3.4460\ 

(4.2) 



B. Intermediate values of M /3 

For > 0.5, the system evolves with damped oscil- 
lations. The transition between non-oscillating and os- 
cillatory mid-late time behavior can be seen in Fig. [5J 
where we show the time evolution of the 4* and Q com- 
ponents with i^2foT 13^ 5-10-3, 0.25, 1 . The behavior 
for higher values of £ is qualitatively similar. 

In this oscillatory regime we have found, for each value 
of M'^P, two modes. In Table U we present the corre- 
sponding QNM frequencies (for £ = 2), computed using 



FIG. 1: (Color online) function of I for different val- 

ues of M^P < 0.4. The dots indicate data from our numerical 
methods, continuous lines indicate the fit (14.211. 
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FIG. 2: (Color online) Time-profiles for the I*!'! (upper panel) 
and |0| (lower panel) components with 1 = 2 and M^fi = 
5-10-^0.25,1. 



both numerical methods described above; we find that 
the agreement between the two approaches is always bet- 
ter then 0.4%. As we discuss in Section flVCi these two 
modes belong to two different branches, which we term 
"gravitational" and "scalar" ; thus we can consider them 
as the "fundamental" modes, i.e. the lowest lying modes 
of these two branches. We stress that these names refer 
to the large /? limit of the modes, but both perturbations, 
^E" and 8, oscillate with both modes^. 

The three different i = 2 modes are shown, for 10^^ < 



^ This happens for M'*/? < 100; for larger values of M'^P, the scalar 
perturbation © oscillates with one mode only, as discussed in 
Section [TVCl 
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TABLE I: Quasinormal frequencies for the oscillatory modes 
with several values of and i = 2. We compare the results 
obtained with the time domain (TD) approach with those 
obtained with the frequency domain (FD) approach. 



TD 


FD 


0.50 0.276 - 0.0967 i 
1.98 - 0.145 i 


0.276 -0.0936 i 
1.97 - 0.144 i 


1.00 0.291 - 0.0970 i 
1.43 - 0.142 i 


0.292 -0.0971 i 
1.43 - 0.142 i 


10.0 0.340 - 0.0980 i 
0.634 - 0.110 i 


0.340 - 0.0983 i 
0.634 - 0.110 i 


100 0.366 - 0.0921 i 
0.501 - 0.0952 i 


0.367 - 0.0919 i 
0.501 - 0.0954 i 


oo 0.374 - 0.0890 i 
0.484 - 0.0967 i 


0.374 - 0.0890 i 
0.484 - 0.0967 i 



M^P < 10^, in Fig. [31 where the dotted-dashed line refers 
to the non-oscillating mode, the continuous line to the 
"gravitational" oscillating mode, and the dashed line to 
the "scalar" oscillating mode. We can see that, for small 
values of /3, the non-oscillating mode Wno, which domi- 
nates the time profile, is excited together with the grav- 
itational oscillating mode; for /3 = 0.3 all three modes 
are present, and for larger values of f3 the two oscillat- 
ing modes are present. Qualitatively similar plots can be 
found for also for £ — 3 and £ = A. The time evolution 
of ^' for = 0.3, which is a combination of the three 
modes, is shown in in Fig. 21 



2.0 - 
S 1.5 - 
^ 1.0 - 



Gravitational mode 
Scalar mode 





10"" 
1 1 III 1 


10-' 


10° 


lo' 

M''p 


10" 


lo' 


10^ 


10 




Gravitational mode 

Scalar mode 

Non oscillatory mode 



lo' 



lo" 



FIG. 3: (Color online) Real (upper panel) and imaginary 
(lower panel) parts of the fundamental QNMs as functions 
of pioT £ = 2. 
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Oscillatory modes: 






Mm ^0.265-0.0911 i 

gov 






Mm ^2.49- 0.172 i 

1 " 1 , 1 


1 , 1 , 1 , i 



150 200 
t/M 



250 300 



FIG. 4: (Color online) Time profile for the component from 
the time-evolution approach (solid line) for M*/3 — 0.3 and 
£ — 2, compared with a combination of oscillatory and non- 
oscillatory modes (dashed line). 



to the fact that the convergence of the iteration scheme 
in the frequency domain approach is more difficult for 
small values of /3. On the other hand, the time-profiles 
are usually available for all the /3 range considered, but 
the extraction of the frequencies from them is not always 
possible. However, we remark that the concordance of 
the two methods is very good in a wide range of param- 
eter space. 



C. Large M*/3 limit 

A time-profile for the wave function for AI'^/S ~ 100 
and £ = 2 is presented in Fig. O The data for the ^ 



9- lO"'' 




10 I — < — r 





10- 1 lo-'^ r I I . I I I I I . 

100 200 300 400 500 100 200 300 400 500 



FIG. 5: (Color online) Time evolution of |^| and |0|, for 
M^P = 100 and £^2. 

component are consistent with a two-mode fit. The val- 
ues obtained are 



It should be mentioned that the numerical determina- 
tion of the different modes for the same value of is 
not an easy task. For instance, neither of the two ap- 
proaches is able to find the scalar non-oscillating mode 
for ^ 0.5. The numerical difficulties are related 



Mwgrav = 0.3736- 0.088991 



0.4837- 0.096711, 



(4.3) 
(4.4) 



which coincide, up to numerical precision, with the 
complex frequencies of the (lowest lying) QNMs of 
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Schwarzschild BHs in Einstein's theory for gravitational 
(wgrav) and scalar (wgc) perturbations [22j. The data 
for the Q component, instead, is consistent with a single 
mode fit, with frequency uj^c- The obtained frequencies fit 
the numerical data very accurately. We can conclude that 
in the /3 — > oo limit and for low multipole numbers £, the 
gravitational perturbations and the scalar field oscillate 
with the QNMs of Schwarzschild BHs: the former, with a 
combination of the scalar QNM and of the gravitational 
QNM; the latter, with the scalar QNM. This behavior 
can be easily understood if we consider the /3 — oo limit 
of the perturbation equations: 



dr'i 



e 



/ 



e{i + 1) 2M 



* = 



967r Mf 



e. 



e = 0. 



(4.5) 
(4.6) 



These equations show that, as discussed in Section |TT1 
the limit /3 — > oo does not correspond to the General 
Relativity limit. Indeed, the gravitational field is coupled 
with the scalar field: equation (|4.5p for ^E* is sourced by 
Q. To recover General Relativity, one should restrict to 
the solutions with = 0; note that 9 = is solution 
of the /? — >■ oo equations (|4.5p . (|4.6p . not of the general 
equations (g^, (j^ . 

Eq. (|4.6p coincides with the equation for scalar field 
perturbations of a Schwarzschild BH in General Relativ- 
ity. It does not depend on 4*, and can be solved sepa- 
rately, yielding the well known scalar QNM frequencies 
of Schwarzschild BHs Once Eq. (j4.6|) is solved, 

one can solve Eq. (|4.5p . treating it like the equation 
of a forced oscillator, since Q{r) can be considered as 
"known" . The homogeneous equation associated to (14.51) 
yields the gravitational QNM frequencies, like Wgrav [221 > 
whereas the source oscillates with frequency Wsc ■ Its solu- 
tion ^(r). at very late times, oscillates with Wsc only, but 
at earlier times it is a combination of the two frequencies, 
as we have found in our numerical integrations. Further- 
more, if O = 0, Eq. (|4.5|) is trivially satisfied, whereas 
Eq. (|4.6p simply becomes the Regge- Wheeler equation 
for gravitational perturbations of a Schwarzschild BH. 
This explains why in the /3 — )■ oo limit both scalar and 
gravitational QNMs are eigenfrequencics of perturbation 
equations. Therefore, no matter how large the coupling 
constant /? is, DGS gravity leaves a peculiar signature in 
the gravitational spectrum of a Schwarzschild BH. The 
actual detectability of this signature is discussed in Sec- 
tion El 



D. Late time power-law tails 

Our results clearly indicate that, for large enough val- 
ues of the coupling constant /3, there is a power-law tail 
dominating the signal of the first multipolar numbers at 



very late times (after the ringdown). Typical time pro- 
files are shown in Fig. [51 The observed late time power- 
law tails arc consistent with the expression 

1>(t,rfi.cd)=t-(^^+^'(;) . (4.7) 

The result (|4.7p can be analytically considered in the 
large r limit. In this limit, the equations decouple and 
previous results in the literature [3l[ are applicable. The 
tails are universal, in sense that they show no depen- 
dence on the parameters M and /3. Note that the same 




FIG. 6: (Color online) Time evolution of I*!'! and |0| for sev- 
eral values of £. Straight lines indicate power-law decay. 

behavior (|4.7p characterizes also Schwarzschild BHs [1^ , 
implying that a gravitational- wave observation of the tail 
would not be able to discriminate DCS gravity from Gen- 
eral Relativity. 



E. Inclusion of a mass term in the Lagrangian 

A relevant question is how the inclusion of a non- 
vanishing potential V{-d) in the action (|2.ip affects re- 
sults discussed above. Here we consider the simplest 
potential, by including a mass term for the scalar field, 
i.e. V{'d) = m^t?^, with m = GM/{hc) and M the physi- 
cal mass of the field. We note that, if we consider a solar 
mass BH, M = 10-^^{mM) MeV; therefore, mM = 
for a massless field, 10^^'^ < mM < 1 for ultra-light ax- 
ions m, mM - lO^s for a pion field, and mM - lO^^ 
for a scalar field at the electrowcak scale. 

We note that the inclusion of a mass term destroys 
the shift symmetry of DCS gravity, i.e., invariance under 
1? — >■ iD + k, with k a constant. If one takes this as a fun- 
damental symmetry, which could presumably be broken 
only at the electroweak scale, it would imply that m is 
of the electroweak size [s^l- Nevertheless, for generality 
we do not impose any a priori constraint on the mass of 
the scalar field. 
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The mass term affects only the perturbation equation 
for the scalar field. In particular only V22 in Eq. p.6p is 
affected, and its general form for massive scalars is 



V2, 



f 



£{£+!) 



1 



5767rA/2 

r6/3 



2M 



(4.8) 

Note that any scalar potential V{i}) whose expansion for 
1? <€. starts at least quadratically, i.e. ¥{19) ^ ■m?-d'^ + 
0(1?^), would give the same potential V22 as Eq. (|4.8|1 . 
This is also the case of the periodic potential V{'d) ^ 
cos d for ultra- light axions [33 . 

Moreover the inclusion of a mass term affects the 
boundary conditions p.lSp for the scalar field. In fact 
at infinity we have 



e 



(4.9) 



Our numerical methods arc capable of computing QNMs 
for massive scalar perturbations whose mass is mM < 
0.2, which includes the case of ultra-light axions. We 
report that numerical results perfectly agree with our 
analytical expectations in Section IIV CI In fact, in the 
large /3 limit, the inclusion of the potential only af- 
fects Eq. (j4.6[) and in turn it modifies only the scalar 
branch of modes: the QNM spectrum consists in the 
usual gravitational modes plus massive scalar modes of 
a Schwarzschild BH. 

For smaller values of the coupling constant {M^(3 < 
100) the analytical limit discussed in Section HV CI breaks 
down and both gravitational and scalar modes are af- 
fected by the scalar potential. Qualitatively, the spec- 
trum for massive scalar perturbations is analogous to the 
one shown in Fig. [5] However, for gravitational modes, 
the dependence on the scalar mass is very mild. The real 
part is almost insensitive to m (at least for mM < 0.2), 
whereas the imaginary part changes as much as 5% for 
M*/3 - 1 and mM - 0.2. Thus, as expected, DCS 
gravity leaves a signature in the QNM spectrum of a 
Schwarzschild BH even if a scalar potential of the form 
(jl.Sp is included. Presumably similar results hold for 
larger values of mM and for more general potentials 



V. DISCRIMINATING THE QNMS: NO-HAIR 
TESTS 

Let us now consider what kind of information one 
can extract from gravitational- wave observations of black 
hole ringdowns, i.e., from the observation of the quasinor- 
mal modes of black holes [H, [s^, HH . 

What we ideally would like to do is to use gravitational- 
wave measurements to test General Relativity and/or 
to rule out alternative candidate theories. The detec- 
tion of two modes in General Relativity would probably 
mean these modes are the £ = 2 and € = 3 fundamen- 
tal modes, with frequencies AIuj — 0.37367— 0.08896 i 
and Mu} ^ 0.59944 - 0.09270 i, respectively [12]. On 



the other hand, two-mode measurements in DCS gravity 
could stand for the lowest £ = 2 modes, which in DCS 
gravity with large M'^P are Mwgrav = 0.3736 - 0.08899 i 
and Mujsc = 0.4837 - 0.09671 i. The question we now 
address is the following: what minimum signal-to-noise 
ratio is required in order to be able to discriminate two 
ringdown signals, and then to test DCS gravity? In other 
words, how can we tell if there really are two or more 
modes in the signal, and can we resolve their parameters? 
If the noise is large and the amplitude of the weaker sig- 
nal is very low, or the two signals have almost identical 
frequencies, the two modes could be difficult to resolve. 
If we can resolve the two modes, then tests of Chern- 
Simons predictions can be performed. 

A crude lower limit on the SNR required to resolve fre- 
quencies and damping times was presented in [22>,34.,35i]. 
The analysis uses the statistical uncertainty in the de- 
termination of each frequency and damping time, which 
a standard Fisher Matrix calculation estimates to be 

illlllll, 



0.1 

par < 65M. 



(5.1) 
(5.2) 



Here, p is the signal-to-noise ratio (SNR), / = Rc[ci;]/27r 
and T = 1/Im[ci;] and ak is the rms error for variable k. 
The numbers above assume white-noise for the detector, 
and equal amplitudes for the two modes. A natural cri- 
terion ( d la Raylcigh) to resolve frequencies and damping 
times is 

1/1 - ./2I > max(cr/i,CT/J , |ti - T2\ > max(crri,crr2) ■ 

(5.3) 

In intcrferomctry this would mean that two objects 
are (barely) resolvable if "the maximum of the diffrac- 
tion pattern of object 1 is located at the minimum of 
the diffraction pattern of object 2". We can introduce 
two "critical" SNRs required to resolve frequencies and 
damping times, 

/ ^ max(pg/^ , pg/J ^ ^ max(pcr^^ , po-^J 

^crit If r \ ' Pcrit i i 

IJI-J2I \Ti-T2\ 

We find the following estimates, 

pLt ~ 6, 
Pent - 150. 



(5.4) 

(5.5) 
(5.6) 



Thus, for SNRs larger than 6, one can distinguish the 
two vibration frequencies in the signal, and is also able 
to discriminate between the General Relativistic and the 
DCS prediction. For SNRs larger than 150, one can also 
measure and discriminate the two different lifetimes. In 
other words, SNRs larger than 6 allow one to discrimi- 
nate between the £ = 3 ringing frequency and the "scalar- 
field-type" gravitational mode in CS gravity with large 
/?. SNRs larger than 150 would allow one to disentan- 
gle even the lifetime of each mode. We also note from 
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Table U and from Fig. [3] that for smaller values of /3 the 
frequency of the (fundamental, "scalar") mode is larger 
and then closer to the £ ~ 3 mode of General Relativity; 
to discriminate between them, a larger SNR would be 
required. 

The results and discussion above assume that both 
modes have the same amplitude. In that sense, the re- 
sults above represent a lower limit for the two modes 
to be discernible. In general the relative amplitude of 
the two modes depends on the physical process exciting 
them and on the coupling parameters of the theory. For 




t/M t/M 

FIG. 7: Dependence of the gravitational-wave signal on the 
relative amplitude of the initial gaussian profile, for /3 — 100 
and £ = 2. 

instance, the relative amplitude has a strong dependence 
on the initial amplitude of each wavepacket, as defined 
by equation (|3J2)) . This is depicted in Fig. [7] for £ = 2 
and f3 = 100. This plot shows that when A1/A2 = 1 
the scalar and gravitational modes compete and the re- 
sult is a damped beating pattern. When A1/A2 = 1000 
the gravitational mode dominates the intermediate-time 
evolution. It would be very interesting to determine the 
relative amplitudes of these modes for physically inter- 
esting situations, like extreme-mass-ratio inspirals, but 
this is outside the scope of the present work. 

VI. CONCLUSIONS 

We have found that Schwarzschild BHs in DCS mod- 
ified gravity arc stable against axial and polar pertur- 
bations. Indeed, an extensive investigation of BH os- 
cillations, performed using two different numerical ap- 
proaches, only yields stable modes, either oscillating or 
non-oscillating. 

Polar perturbations obey exactly the same master 
equation as in General Relativity, and therefore BHs in 
DCS gravity oscillate at the same polar frequencies. Ax- 
ial perturbations, instead, couple to a scalar field, en- 
larging the spectrum of resonances in the gravitational 



sector. In particular, the ringdown of a BH in DCS grav- 
ity is a superposition of two different QNM sectors. For 
large values of the constant /3, which is associated to the 
dynamical coupling of the scalar field, one of these sectors 
corresponds to the gravitational and the other sector to 
scalar-field QNMs of Schwarzschild BHs in General Rel- 
ativity. Thus, a golden opportunity to test these theories 
is by detection of BH ringdowns. As shown in Section fVl 
a modest SNR (> 6) could be sufficient to discriminate 
between General Relativity and DCS modified gravity. 
These estimates assume very special relative amplitudes 
between the modes. Accurate estimates, as well as con- 
straints on the coupling parameters, require the calcu- 
lation of accurate waveforms for physically interesting 
processes exciting these ringdown modes. 

The problem dealt with here is also interesting for a 
number of other reasons, in particular because we expect 
such kind of problems, i.e. QNMs described by a system 
of coupled second order DDEs, to be a general feature of 
alternative and more intricate theories; surprisingly there 
are very few studies of this kind of system in General 
Relativity. 

Finally, we detail in Appendix |X] how ghost- 
instabilities develop in this theory when /3 < 0, by a 
careful analysis of the instability timescale and other fea- 
tures. 

Generalization of our results to rotating black holes 
is of utmost importance, given that many astrophysical 
black holes are rapidly rotating. Rotating solutions in 
DCS gravity are only partially understood [l^ [s!], we 
hope to come back to this issue in the near future. 
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Appendix A: Ghost-like instabilities for /3 < 

In this Appendix we study unstable modes for the sys- 
tem (|2.8p . (|2.9p . which arise for (3 < 0, i.e. for the wrong 
sign for the kinetic energy in the action (j2.ip . In par- 
ticular we discuss a peculiar instability, arising at large 
multipoles £, which we believe may be seen as a general 
signature of ghost-like instabilities at linear level. 

For /? < 0, our numerical approaches both show that 
the amplitudes of the gravitational and scalar field grow 
exponentially with time: the spacetime is unstable. The 
agreement between the two methods is excellent (to 
within the last significant digit), thus results presented 
here can be reproduced by both methods. 

For small values of Af*^!/?! the growth is purely expo- 
nential, ~ e"°°*. The non-oscillatory exponential coeffi- 
cient Wno depends on /3 and i, as presented in Table HIl 



Expression (jA3[) implies that the complete perturbation 
(taking into account all multipole components) is always 





TABLE II: Non-oscillatory exponential coefficient Mojn 
several values of and £. 





£ = 2 




^ = 3 




M*p = -1 








Mujno 


I 




-0.05 


5.894 


-0.05 


8.391 


2 


1.115 


-0.10 


4.111 


-0.10 


5.871 


3 


1.629 


-0.50 


1.706 


-0.50 


2.467 


4 


2.142 


-1.00 


1.115 


-1.00 


1.629 


5 


2.655 


-2.00 


0.6666 


-2.00 


0.9930 


10 


5.215 


-3.00 


0.4382 


-3.00 


0.6710 


20 


10.31 


-4.00 


0.2650 


-4.00 


0.4358 


30 


15.30 


-4.791 


0.0547 


-5.233 


0.0752 


50 


24.69 



FIG. 8: (Color online) Evolution in time of |*[ and |9| for 
for 1 = 2 and negative values of M'^p. 



unstable: for any /3 < there is always a multipole £ for 
which oj/ > 0. 

In Tabic IIIII we present the fundamental, unstable 
mode, for large values of (3 and different values of £. The 
imaginary part of the fundamental unstable mode grows 
linearly with £, i.e. the instability timescale decays lin- 
early with £. 



MUJy, 



-0.22- 0.19£- 



0.30 + 0.69^ 

(Af4 1^1)0-45 



(Al) 



for any £ and small enough M4|/3|. We notice that the 
instability timescale t = 1/Im[w] is shorter (the instabil- 
ity is stronger) for smaller /3 and for larger £. From the 
expression above one expects that, for large enough |/3|, 
pure exponentially-growing modes cease to exist (uJno is 
negative for large enough |/3|). In fact if /3 < — |/3no|, the 
late time decay is dominated by an oscillatory exponen- 
tial mode. This is depicted in Fig. [8] 

Furthermore our results show that, for fixed £ and as /3 
further increases, there exists a critical value /3crit, such 
that 



TABLE III: Fundamental unstable mode for different values 
of P and taking into account the multipole components up to 
^ = 30 



Fundamental unstable mode, 



■ UJR + IUJI 





M-^P 


= -10 


M^P 


= -20 


M^P 


= -30 


M^P -- 


= -40 


£ 


MojR 


Muji 


MujR 


Muji 


MujR 


Muji 


Mu}R 


Muji 


3 


0.5387 


0.0034 


0.5835 


0.0060 










4 


0.7342 


0.0540 


0.7771 


0.0371 


0.7973 


0.0177 


0.8074 


0.0037 


5 


0.9154 


0.0964 


0.9656 


0.0676 


0.9885 


0.0443 


1.0002 


0.0278 


10 


1.7929 


0.2884 


1.8862 


0.2186 


1.9251 


0.1732 


1.9456 


0.1422 


20 


3.5266 


0.6625 


3.7061 


0.5188 


3.7794 


0.4291 


3.8181 


0.3686 


30 


5.2563 


1.0362 


5.5216 


0.8186 


5.6295 


0.6848 


5.6865 


0.5947 



Im[w(/3crit)] =0, (A2) 

and for (3 < — |/3crit| the modes change from unstable to 
stable. This critical value depends on £ and its depen- 
dence is very well fitted by a quadratic function 

= -2.77^2 _ (A3) 



Generically our results imply that for any value of 
/? < there is an instantaneous instability which de- 
velops once all the multipolar components are taken into 
account. This is related to the choice of the wrong sign 
for the kinetic term of the scalar field in the action, and 
it is the signature of ghost-like states at the linear level. 
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